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Abstract 

Background: In Japan, a shortage of physicians, who serve a key role in healthcare provision, has been pointed out 
as a nnajor nnedical issue. The healthcare workforce policy planner should consider future dynamic changes in 
physician numbers. The purpose of this study was to propose a physician supply forecasting methodology by 
applying system dynamics modeling to estimate future absolute and relative numbers of physicians. 

Method: We constructed a forecasting model using a system dynamics approach. Forecasting the number of 
physician was performed for all clinical physician and OB/GYN specialists. Moreover, we conducted evaluation of 
sufficiency for the number of physicians and sensitivity analysis. 

Result & conclusion: As a result, it was forecast that the number of physicians would increase during 2008-2030 
and the shortage would resolve at 2026 for all clinical physicians. However, the shortage would not resolve for the 
period covered. This suggests a need for measures for reconsidering the allocation system of new entry physicians 
to resolve maldistribution between medical departments, in addition, for increasing the overall number of clinical 
physicians. 

Keywords: Forecasting the number of physicians. Shortage of physician. Maldistribution, System dynamics 



Introduction 

Health policy planning for human resources has become 
an internationally high priority [1]. In Japan, a shortage 
of physicians, who serve a key role in healthcare provision, 
has been pointed out as a major medical issue [2]. The 
Ministry of Health, Labour and Welfare (MHLW) 
reported that there is lack of absolute and relative num- 
bers of physicians, as well as maldistribution among re- 
gions and medical departments. To ease the personnel 
shortage in Japan, the MHLW estimated that hospitals 
need a further 24 000 physicians [3]. 

It has been suggested that the policy of decreasing 
medical school enrolment quotas is a major factor be- 
hind the absolute shortage of physicians [4,5]. In 1982 
and 1999, the policy was approved by the Japanese 
Government to prevent future oversupply. The medical 
school enrolment quota had decreased from 8 280 
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students in 1984 to 7 625 in 2007. However, the number 
of physicians per unit of population in Japan is low com- 
pared with other developed countries [6]. According to 
Organization for Economic Co-operation and Develop- 
ment (OECD) data, the number of medical doctors per 1 
000 persons in Japan was 2.2 in 2008; this number was 2.4 
in the US in 2009 and 2.7 in the UK in 2010. The number 
of physicians in Japan has always been lower than the 
OECD average. In 2008 the MHLW identified the phys- 
ician supply issue, and recognized that there is not only 
maldistribution, but also a shortage of absolute numbers 
of physicians [7]. To address the shortage, the Japanese 
Government decided to immediately increase the medical 
student enrolment quota. The Democratic Party of Japan 
(DPJ), which was elected into government from 2009 to 
2012, promised to increase the quota to 1.5 times what it 
was in 2009 in order to raise the number of physicians per 
unit of population to the OECD average level [8]. 

Maldistribution of physicians at both the regional and de- 
partmental level has been reported by the MHLW [3]. Re- 
gional maldistribution is identified as a factor contributing 
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to the generation of areas that fail to meet the basic 
healthcare needs of their community. Maldistribution 
among medical departments causes insufficient provision 
of some medical care. In Japan, obstetrics/gynaecology 
(OB/GYN) and emergency departments are especially 
mentioned as being in need of more physicians than other 
departments. Specifically, the number of OB/GYN physi- 
cians is decreasing, and the department is facing future 
shortages [9]. In Japan, medical school graduates determine 
their own speciality and location, so physician maldistribu- 
tion may worsen. To address regional maldistribution, a 
bonding scheme has been proposed to ensure that gradu- 
ates work in specified regions. No comprehensive measures 
have been proposed to address medical department-level 
maldistribution, and it is at risk of worsening in the future. 

The healthcare workforce policy planner should con- 
sider future changes in physician numbers. Effective 
planning for the future medical workforce has become in- 
creasingly necessary in Japan because of changes in and 
uncertainty of physician supply in the face of the popula- 
tions requirement for their services [10]. Accurate fore- 
casting of the workforce supply is essential for effective 
planning. Forecasting methods include the need-based 
model, the demand-based model and benchmarking 
[11-13]. These methods are limited by assumptions about 
changes in physician availability across incomparable 
communities or areas [11]. We argue that the forecasting 
approach used should be able to describe the dynamic 
variables involved in the physician supply system and have 
utility even if the needs of the community and health 
system differ. 

The system dynamics approach permits the inclusion 
of a dynamic factor in the models. We modeled the 
physician supply system to estimate the future number 
and sufficiency level of clinical physicians and forecast 
the absolute number of physicians required. Furthermore, 
we compared the sufficiency levels between all clinical 
physicians and specialists to evaluate department-level 
maldistribution. The purpose of this study was to propose 
physician-supply forecasting methodology by applying SD 



modeling to estimate future absolute and relative numbers 
of physicians. 

Materials and methods 

Data collection 

We obtained data on the number of physicians from the 
national surveys of physicians, dentists and pharmacists 
conducted from 1972 to 2008 by the MHLW [9]. We 
extracted data on physician numbers, age, gender, speci- 
ality, and type of workplace from the survey. Medical 
specialities were defined according to the MHLW classi- 
fication. To account for increases in medical school en- 
rolment, we obtained data from the Ministry of 
Education, Culture, Sports, Science and Technology 
(MEXT) [14]. In addition, we used data from the basic 
school survey carried out by MEXT to identify past 
graduation numbers, and trends in speciality choice after 
graduation [15]. Moreover, the retirement rate was cal- 
culated by trend data considering the age-sex pyramids 
of the medical population. Other data sources are sum- 
marized in Table 1. 

SD modeling approach 

Our model was based on SD, a methodology and com- 
puter simulation modeling technique for framing, under- 
standing and discussing complex issues and problems. 
The approach was developed by computer pioneer Jay 
W Forrester in the mid-1950s [16]. It addresses the flow 
of people, processes, materials and information by ex- 
posing a tenet of all complex systems - feedback loops. 
According to SD, these loops are the major drivers of 
system behaviour. This methodology also enables us to 
address non-linear SD, which governs many real-life 
phenomena [16]. The approach has been used in a 
variety of contexts, including human resources planning, 
to gain an understanding of a system with complex dy- 
namic and nonlinear interacting variables [17,18]. In- 
ternal feedback systems are represented in the SD 
model. The model is an interlocking set of differential 



Table 1 Data sources for the simulation model 



Variables 



Data sources 



Current workforce in baseline year (2008) 
Medical school quota (current) 
Pass rate at the national examination 
Career option rate 
Retirement rates 

(include: death, change of speciality) 
Population projection in Japan 

Selection speciality rates 



The number of physician reported by the Ministry of Health, Labour and Welfare in 2008 
Bulletin about scheme of increases medical school quota 
Announcement about national examination for medical doctors 

The school basic survey carried out by the Ministry of Education, Culture, Sports, & Technology 
We set the constant rate from trend data for the numbers of physicans reported by the MHLW 

Population projection for Japan 

(National Institute of Population and Social Security Research) 
Questionnaire survey about clinical training in 2010 
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Figure 1 Conceptual scheme of the Japanese physician career path. 



algebraic equations developed from a broad spectrum of 
relevant measured and experiential data [16]. 

We used specialized software, STELLA® ver. 8.1.1 (isee 
systems, Lebanon NH, USA) to implement the SD 
models [19]. 

Causal loop and stock and flow 

A causal loop diagram identifies the structures and inter- 
action of feedback loops, and consists of variables for 
causal links. A causal link connects a cause variable with 
an effect variable by a link with a positive or negative 
change. To construct the physician supply model, we based 
the causal loop on the Japanese physician career path 
programme (Figure 1). 

System accumulation is described by stock (or level) 
variables, such as the number of physicians, while the 
rate of change is described by flow variables, such as the 
annual number of physicians who retire. The net inflow 
determines the rate of stock change, that is, its time de- 
rivative. Therefore, the interaction between stock and 
flow is calculated mathematically by the following differ- 
ential equation. In the equation, inflow and outflow de- 
note the values of the inflow and outflow at any time 
between the initial time and the present time t [19]. 

^^^^^^^^ = Inflow{t)-Outflow{t) (1) 



Uncertainty projection 

Uncertainty in health projections must be assessed so 
that planners can anticipate possible variations and 
adapt the planning of human resources in consequence 
[10]. We conducted a sensitivity analysis to estimate the 
influence of the uncertainties of medical school quota 
and trends in graduate speciality choice on the forecast 
sufficiency level. We set the following quota and special- 
ity choice scenarios from previous variables: increasing 
or decreasing the quota from the current number by 
15%, and increasing or decreasing the number of gradu- 
ates from the current number choosing OB/GYN by 1%. 

Criterion for evaluation the number of physicians 

In planning the health workforce, a supply-need gap 
analysis must be conducted [11]. In Japan, there is no 
criterion for evaluating whether the number of physi- 
cians is sufficient or not. Conventionally, the number 
has been evaluated by a comparison with the OECD 
average level. In this evaluation, health workforce plan- 
ners cannot account for the number of physicians that 
are required by patients or medical workers. We defined 
the physician number sufficiency level using data from 
the national survey carried out by the MHLW in 2010 
to enable us to evaluate whether the future physician 
numbers are sufficient or not. 
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Table 2 Key assumption of the simulation model 

Variable Key assumptions 



New medical graduates 



Retirement 



Domestic students only 

The number of graduates is equal to the 
medical student quota 

Retirement, deaths and movement speciality 
rates remain constant 



sufficiency level 



the forecasted number of physician 
the required number of physician 
X corrective coefficient 

(2) 



In equation (2), the required physician number is de- 
rived from the claims of hospital administrators, and 
does not include clinics. The corrective coefficient cor- 
rects for missing clinic data and is the ratio of physicians 
who work in a hospital to those who work in a clinic. 



In this study, if the sufficiency level is >1, we consider 
that the physician number is sufficient. When the suffi- 
ciency level is <1, we consider that there is a physician 
shortage. Furthermore, we evaluated speciality-level mal- 
distribution by the difference in sufficiency level between 
all clinical physicians and OB/GYN physicians. If the 
sufficiency level of OB/GYN is <1 and the sufficiency 
level of all clinical physicians is >1 (and vice versa), we 
evaluated that there is speciality-level maldistribution 
between OB/GYN and other specialities. 

We constructed two models using the SD approach: 
(1) total number of clinical physicians, and (2) number 
of OB/GYN physicians. The base year was 2008 and the 
simulation is forecast up until 2030. The outcome 
measures were the numbers of clinical physicians and 
specialists (head count and per 1 000 persons) and suffi- 
ciency level. In order to test the validity of the models, we 
conducted a retrospective comparison of forecast and ob- 
served numbers in 1986 to 2008. The SD model method- 
ology used for the health workforce forecasting in this 
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study is displayed in Figure 2. The key model assumptions 
are summarized in Table 2. 

The validity of the SD model 

The validity of model forecasts should always be tested. 
As a necessary step in SD methodology to validate the 
model, historical and simulated data should be compared 
for specific years [20]. We examined the simulated and 
observed physician numbers in 1998, 2000, 2002, 2004, 
2006, and 2008. The method is to compare the simulation 
data and observed historic data, and calculate the relative 
and mean squared errors of the number of physicians: 



\iyt-yt)\ 
yt 



MSE 



(3) 



(4) 



Here, yt represents the observed number of physicians 
in the year t; represents the simulated number of phy- 
sicians in the year e is the relative error of the number 
of physicians and mean squared error (MSE) is the mean 
of the squares of the relative errors. 

Results 

Figure 3 displays the causal loop and stock and flow dia- 
grams of the SD model. There are six main sectors in 
this diagram; enrolment, national examination, clinical 
training, clinical practice, choice of speciality, and retire- 
ment. Figures 4 and 5 display the forecast numbers of all 
clinical physicians and OB/GYN specialists. In 2008, 



there were 271 897 clinical physicians in Japans medical 
workforce. Of these, 10 388 were OB/GYN specialists. 
The number of clinical physicians was forecast to in- 
crease by 370 345 during 2008 to 2030. This represents 
an average annual growth rate of 1.6%. Figure 6 displays 
the forecast number of physicians per 1 000 persons and 
the population projection for Japan by the National 
Institute of Population and Social Security Research in 
2012. The number of physicians per 1 000 persons was 
2.2 in 2008, and was forecast to increase by 3.18 per 10 
00 persons during 2008 to 2030. This represents an aver- 
age annual growth rate of 2.2%. Figure 7 displays the 
forecast sufficiency level. The number of clinical physi- 
cians was in shortage in 2008. However, it was forecast 
to gradually increase and reach a sufficient level by 
2026. The number of OB/GYN specialists was forecast to 
increase by 13 498 by 2030. The forecast sufficiency level 
of OB/GYN specialists was less than one during the entire 
period 2008 to 2030. Therefore, we evaluated that there 
would be speciality-level maldistribution. The sufficiency 
level growth rate for all clinical physicians was 36.21% 
during the forecast period. On the other hand, the rate for 
OB/GYN specialists, who are especially required, was only 
30.1%. The sufficiency level of all clinical physicians sur- 
pass one, while the sufficiency level of OB/GYN specialists 
would not have reached one after 2026. 

The results of the sensitivity analysis are shown in 
Figures 8 and 9. The analysis of all clinical physicians 
under the increased quota scenario shows an accelerated 
growth of sufficiency levels from shortage to sufficient. 
In the decreasing quota scenario, the growth of suffi- 
ciency levels was delayed compared with the baseline 
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Figure 4 Forecast numbers of all clinical physicians, 2008 to 2030. 
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Figure 5 Forecast numbers of obstetrics/gynaecology specialists, 2008 to 2030. OB/GYN, obstetrics/gynaecology. 



scenario, and did not reach one. For OB/GYN special- 
ists, neither the increased nor decreased quota scenario 
leads to a sufficient level (Figure 8). When comparing 
the influence of the medical quota versus the speciality 
option, the latter had a more significant impact on the 
sufficiency level than the former (Figure 9). 

Simulated and observed historical data for specific 
years were compared to validate the SD model. Table 3 
presents this comparison as well as the calculated MSE. 



Discussion 

Our model indicates that the number of clinical physi- 
cians in Japan will continue to grow during 2008 to 
2030. However, under the baseline scenario, a sufficient 
number of physicians would not be reached until 2026. 
The DPJ promised in its manifesto to raise the number 
of physicians per unit of population to the OECD 
average level. We determined that the number of clinical 
physicians per 1 000 persons would not surpass the 



3.5 



1.5 



1.0 



I 0.5 



0.0 



-The number of physician per 1000 
-population 




YEAR 

Figure 6 Forecast numbers of all clinical physicians per 1 000 persons. 



80,000 "5 

a 
o 

Q. 

60,000 in 

c 



.o-v^ J> ^<S^ .o^^ J' J" 



Ishikawa et al. Human Resources for Health 2013, 11:41 
http://www.hunnan-resources-health.conn/content/1 1/1/41 



Page 7 of 10 



1.5 
1.4 
1.3 



_ 1.2 
> 

^ 1.1 









• • all clinical physicians 

• • OB/GYN 












'^<^^ 



YEAR 



Figure 7 Forecast sufficiency level, 2008 to 2030. 



OECD average until 2029 (the OECD average was 3.1 in 
2008; our forecast number is 3.12 by 2029). Our analysis 
shows that the DPJ promise would be fulfilled within the 
years covered in the SD model. 

Our model forecasts that the shortage in the abso- 
lute number of clinical physicians would be resolved 
during the coverage timespan of the SD model. How- 
ever, the number of OB/GYN specialists would still 



remain in shortage in 2030. Given from above, we 
evaluated that there would be speciality-level maldis- 
tribution. We speculate that the difference in the 
sufficiency-level growth rate raises/expands the maldis- 
tribution. We propose that measures for increasing the 
numbers of OB/GYN specialists is required in addition 
to measures for increasing the overall number of 
clinical physicians. 
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Figure 8 Sensitivity analysis of the sufficiency level of all clinical physicians and obstetrics/gynaecology specialists, 2008 to 2030. 
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Figure 9 Comparison of impact for medical schiool quota and speciality choice on sufficiency level. 



The sensitivity analysis provided an evaluation of the 
impact of each factor. Figure 8 shows that the number 
of physicians would be sufficient in 2023 in the in- 
creased medical student enrolment quota scenario 
(+15%), whereas it would remain in shortage in 2030 in 
the decreased medical student enrolment quota scenario 
(-15%). Figure 9 shows that increasing the number of 
graduates choosing OB/GYN has a greater impact on 
the OB/GYN sufficiency level than increasing the med- 
ical student enrolment quota for OB/GYN. 

The sensitivity analysis, looking from different angles, 
can be interpreted as a scenario analysis. In Japan, most 
physicians build their career as indicated in Figure 1. 
With the matter in mind, the rates of speciality choices 
of potentially causes an increased shortage of the num- 
ber of specialists, as shown in Figure 9 under the -2% 
scenario. In this way, even if the absolute shortage is re- 
solved, the uncertainty of variation in choices of special- 
ity possibly causes shortages in specific speciality, and 
the speciality-level maldistribution. Medical school grad- 
uates voluntarily determine their speciality in Japan. This 
generates the uncertainty of variation in the rates of spe- 
ciality choices. Therefore, we highlight the necessity of 
strategic planning for speciality options, for example, by 
obligating a proportion of medical school enrollees to 
choose a specific department. 

The validity of the model forecasts should always be 
tested. Our model was validated by calculating the relative 
error and MSE from simulated and observed historic data. 
Generally, MSE <10% is considered reasonable [20]. 
According to equations (3) and (4), our MSE was calcu- 
lated as 0.007 <0.1 (Table 3). Furthermore, the relative 



error had no effect on our evaluation of the number of 
physicians. We evaluated the impact of relative error as 
being very low, and therefore concluded that our model 
was valid. 

Some previous studies and the government have also 
estimated the numbers of these professionals by different 
methods. But these estimates have been based on the 
static model method. We performed the estimate consid- 
ering the dynamic change of some value that impacts the 
physician number. As with any modeling exercise, our 
model was constructed under the assumption the system 
would not change. A change in the physician supply sys- 
tem would require further modeling. Therefore, we should 
be rethink the forecasting model structure when the sys- 
tem changes. 

In this study, we defined a criterion for judging whether 
the number of physicians is sufficient or not. The criterion 
was calculated based on the survey performed by the 
MHLW in 2010 [3]. The data were collected using a self- 
report questionnaire for the medical institutions. They 
subjectively provided their required number of physicians 
in 2010, judging from daily clinical practice. Therefore, an 
evaluation using the required number of physicians from 
the survey would be based on the assumption that 



Table 3 The comparison of historic data and simulation 
data 



Year 


1998 


2000 


2002 


2004 


2006 


2008 


MSE 


Historic data 


248.61 


255.79 


262.68 


270.37 


277.93 


286.7 


0.007 


Simulation data 


250.39 


257.88 


264.74 


271.16 


277.09 


283.71 




Relative error 


0.007 


0.008 


0.008 


0.003 


0.003 


0.01 
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answers to the questionnaire reflect regional need. We as- 
sume the required number reflect health statistics, mor- 
bidity rate and prevalence rate and so on, based on their 
honest answers. In the OB/GYN domain, in particular, the 
female-specific influences of fertility rate and the inci- 
dence of gynaecological carcinoma for instance, was a 
consideration in the OB/GYN specialists' opinions. The 
required physician number by MHLW 2010 did not 
provide the demographics of age or gender. Therefore, 
the criterion does not reflect a difference between gen- 
ders. We need to explore the quantification method 
considering the difference. Moreover, the criterion can- 
not take account of the impact of dynamic change on 
need. Generally, as time goes on, it is expected that the 
required number of physicians will change. We need to 
account for this change when estimating the sufficient 
number of physicians needed in the future. 

Our model did not account for the shortage or maldis- 
tribution between regions. SD modeling permits the 
addition of a physician transfer factor to the model. By 
adding this factor, the forecast number of physicians 
would enable us to evaluate the future supply-need 
balance while focusing on the presence or absence of 
regional or departmental shortages. We demonstrated 
that the SD modeling approach can be used to forecast 
changes in the number of physicians and the sufficiency 
level. The model provided two key suggestions. First, the 
absolute shortage of physicians will resolve during the 
forecast period because the number per population 
would reach the 2009 OECD average and the sufficiency 
level would reach one. Second, the number of OB/GYN 
specialists is at risk of remaining in shortage. Measures 
to address the uncertainty in medical trainees' choices 
of speciality are necessary to prevent an increase in 
department-level maldistribution. Obligating a propor- 
tion of medical school enrollees to choose a specific de- 
partment will remove this uncertainty. 

A proactive approach to the planning and manage- 
ment of the medical workforce, involving continued 
monitoring, would facilitate more frequent and less dra- 
matic adjustments to supply. This in turn would reduce 
the likelihood of extreme supply conditions (shortages 
or surpluses) and assist in more effective delivery of 
medical care. SD modeling enabled us to analyze mul- 
tiple perspectives, considering the causal links, factor 
uncertainties, and different scenarios. However, this 
study has some limitations derived from the difficulty in 
obtaining some data. Increasing the coUection of data 
on physician supply and the required number of physi- 
cians would minimize these limitations. 

Conclusion 

In this study, we defined a criterion for judging whether 
the number of physicians is sufficient or not. The 



criterion was calculated based on the survey performed 
by the MHLW in 2010 [3]. Generally, as time goes on, it 
is expected that the required number of physicians wiU 
change. We need to account for this change when esti- 
mating the sufficient number of physicians in the future. 
Our model did not account for the shortage or maldis- 
tribution between regions. SD modeling permits the 
addition of a physician transfer factor to the model. By 
adding this factor, the forecast number of physicians 
would enable us to evaluate the future supply-need 
balance while focusing on the presence or absence of re- 
gional or departmental shortages. We demonstrated 
that the SD modeling approach can be used to forecast 
changes in the number of physicians and the sufficiency 
level. The model provided two key suggestions. First, 
the absolute shortage of physicians will resolve during 
the forecast period because the number per population 
would reach the 2009 OECD average and the sufficiency 
level would reach one. Second, the number of OB/GYN 
specialists is at risk of remaining in shortage. Measures 
to address the uncertainty in medical trainees' choices 
of speciality are necessary to prevent an increase in 
department-level maldistribution. Obligating a propor- 
tion of medical school enroUees to choose a specific 
department will remove this uncertainty. A proactive 
approach to the planning and management of the med- 
ical workforce, involving continued monitoring, would 
facilitate more frequent and less dramatic adjustments 
to supply. This in turn would reduce the likelihood of 
extreme supply conditions (shortage or surplus) and 
assist in more effective delivery of medical care. SD 
modeling enabled us to analyze multiple perspectives, 
considering the causal links, factor uncertainties, and 
different scenarios. However, this study has some limita- 
tions derived from the difficulty in obtaining some data. 
Increasing the coUection of data on physician supply 
and the required number of physicians would minimize 
these limitations. 
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